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We investigate collective excitations coupled with monopole and quadrupole oscillations in two- 
component fermion condensates in deformed traps. The frequencies of monopole and dipole modes 
are calculated using Thomas-Fermi theory and the scaling approximation. When the trap is largely 
deformed, these collective motions are decoupled to the transverse and longitudinal breathing os- 
cillation modes. As the trap approaches becoming spherical, however, they are coupled and show 
complicated behaviors. 
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I. INTRODUCTION 

Since the realization of the Bose-Einstein condensed (BEC) atomic gases , there has been much 
interest in ultracold trapped atomic systems to study quantum many-body phenomena. Besides the 
Bose-Einstein condensates (BEC) , one can now study degenerate atomic Fermi gases 4] and 

bose-fermi mixtures Q. These systems offer great promise to exhibit new and interesting phenomena of 
quantum many-particle physics. 

Sogo and Yabu Q studied the ground state properties of fermi gas composed by two spin components 
with various repulsive interaction and showed the phase transition between the paramagnetic and ferro- 
magnetic states. Furthermore we extended this study on the spin phases of the ground state to that for 
the asymmetric gas with a nets spin in the presence of an external magnetic field and gave the phase 
diagram about these spin phases 

An important diagnostic signal for many-body systems is the spectrum of collective excitations. Such 
oscillations are common to a variety of many-particle systems and are often sensitive to the interaction 



and the structure of the ground state. In ref.Q we have also studied the spin excitation on the dipole and 
monopole oscillations in the two-component fermion condensed system. The frequencies of the collective 
oscillations with the out-of-phase show the dependence on the fermion-fermion interaction while those 
with the in-phase are not sensitive to the interaction. Particularly the frequencies of the out-of-phase 
modes show rapid changes in the phase transition between the paramagnetic and ferromagnetic spin- 
phases. 

The monopole oscillations can be defined only in the spherical symmetric system. In usual experiments 
traps have elliptically deformed shapes, where the monopole oscillations are coupled with the quadrupole 
oscillations. If the trap is largely deformed with the axial symmetry the oscillations must be decoupled 
into the modes of oscillations along the transverse and longitudinal directions; we shall call them the 
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transverse breathing (TB) and longitudinal breathing (LB) modes. 

In this paper we study the coupled oscillation states of the two component fermi gases between the TB 
and LB modes in deformed traps. The oscillation behaviors are expected to show various features as the 
deformation of the trap is changed. 

The underlying theoretical tool to treat dynamic problems of dilute quantum gases is the time- 
dependent mean field theory. This is reduced to the random-phase approximation (RPA) for small 
amplitudes, and the theory in this form has been applied to density oscillations of these systems 
When the single-particle spectrum is regular, the long-wavelength excitations are collective and simpler 
methods can be used to calculate the frequencies, in particular with sum rules lllj or the scaling ap- 
proximation. In this work we calculate the frequencies of these collective oscillations by using the scaling 
method [T3. Ilj. This scaling method is ph ysic ally quite transparent, but it is in fact equivalent to the 
theory based on energy- weighted sum rules 13 1 , as used for example in ref . fll| . 



We will only consider here the case of positive scattering lengths which correspond to a repulsive 
interaction. When the interaction is attractive, the system is superfluid in its ground state and the 
excitation properties are controlled by the energy gap. In addition we will focus only on the collective 
oscillations in the paramagnetic spin phases, and do not deeply discuss those in the ferromagnetic regime. 

In the next section we will explain our formulation, and show calculational results in Sec. 3. In Sec. 4 
we will summarize our work. 

II. FORMALISM 

We begin with the expression for the energy of a dilute trapped condensate in the mean-field approxi- 
mation: 

„ 1 O A 1 — 



s=l,2 j=l 



with 



Tsj ~ L> drj Or, 1 

occ 
n 

where ipns are orbital wave functions indexed by orbit n and component s = 1,2, Ql,t are the longitudinal 
and transverse frequencies of the trapping field and gp is the coupling strength of a contact interaction. 
Here we assume that all fermions have the same mass, m, but we do not need to treat these components 
as the spin state of the same fermion. 

In order to study the collective oscillation, we introduce the following scaling for the fermion wave- 
function: 

^,.(r, r) = e^'V*.. M+^.m^o) ( e ^.Wr r ; e^^r 3 ), (4) 
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where ipn}s is the wave- function in the ground state, \t,l( s ) are the time-dependent collective coordinates 
for the oscillation, and \t,l{s) are the time-derivatives of Xt,l{s) ■ The factor exp(i£) is the Gallilci 
transformation factor which is determined to satisfy the continuum equation. 

Let us obtain the total energy under the above scaling by substituting the wave-function into the 
total energy functional (JJJ. The total energy becomes 

s=l,2 J 

+\™ 2 E / d 3 r{e-^nUrl + rl) + e-^nlri)}p^(r), 

s=l,2 

+ge ^T,i+\L.i+2\T,2+\L,2 J d 3 r P i ° ) (e XT - 1 r T ;e Xl - 1 r3)p i ° ) (e XT - 2 rT;e XL ' 2 r 3 ) (6) 

where the superscript (0) indicates t s j and p s at the ground state. 

In this work we use the Thomas- Fermi (TF) approximation to evaluate the above total energies. The 
TF approximation makes the spherical symmetric in the momentum distribution at each position and 
gives the following relation: 

* ) =Ti5=riS = i l ; (6^)lp?)*. (7) 

In addition the ground state densities ps become functions of flr( r i + r l) + ^L r | ■ 

Here we can then make a change of variables to simplify the appearance of the Thomas- Fermi equations 
as follows. 

* = <^ <i=l~3>, 

- 4m 12 g§tl L n 2 T p 

fJ = f (8) 

where go is a certain coupling constant which is chosen later. 

With the scaled variables and the Thomas-Fermi approximation, the expression for the total energy at 
the ground state is 



I d 3 x { i^nl + x 2 n s ) + ,gnin 2 }, (9) 

« — io 
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where x 2 = \x\ 2 . The TF equations for the densities n\^, are derived by variation of the energy with a 
constrain on the fermion numbers. This yields 

2 

n'( + gn,2 = e\ — x 2 — e/ + B — x 2 

2 _ 

n| + gni = e 2 — x 2 = ej — B — x 2 (10) 



In this equation e\ and e 2 are the Lagrange multiplier to constraint the numbers of the fermion- 1, Ni 
and -2, iV 2 , respectively, and they are rewritten with e/ and B; these e/ and B have the meanings of the 
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Fermi energy and the external magnetic field in fcrmion gas when the two components of fermions are 
described with two different spin states of one fermion. 

In this expression the ground state is determined by the three variables, g, e/ and B, but the two 
variables, g and e/, are not independent because of the scale symmetry in ea. (|10|l HQ. In this work 
we choose go as a critical coupling constant of the paramagnetic and ferromagnetic phase transition at a 
fixed chemical potential in the symmetric system = p 2 . With this choice the fermi energy is fixed 
to be e/ = 20/27. The solution of these equations is discussed in detail in ref. |7j. With fixed values of 
ef = 20/27 and B — 0, the spin phase of the two component fermi gas turns to be 'paramagnetic' (g < 1) 
, 'ferromagnetic \g > 1) in order as the coupling constant, g, increases. 

Furthermore the total energy in collective oscillations is given as 
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with 



f s = - d 3 x ni{x) 



X s = d x x n s (x) 



(12) 



where x 2 = \x\ 2 , lot = Qt/^f and lol = ^l/^f (^f = v^r^^)- The interaction energy Vff appears 
Vff = ffe 2 ^T,i+AL.i+2A T , 2 +2A L , 2 I d^r nx{e XT ' x x T ,e XL < x xz)n 2 {e XT ' 2 XT,e XL *xz) (13) 



In order to describe the collective oscillations, we obtain the variation of the total energy up to the 
order of 0(A 2 ) as 
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with 



= gjd'x (xni^j , (17) 
V 2 = g J (x^n^j , (18) 



V 3 



where the vector A is defined as *A = (At,i, ^t,2, -Vl,2)- The terms propritional to A disappears 
because of the virial theorem. 

Then the classical equation of motion for A is harmonic, giving rise to the following eigenvalue equation 
for the oscillation frequencies, 

[Blu 2 - C] A = 0. (20) 

Before performing actual calculations, we should discuss the collective oscillations in some extreme 
conditions. 

In the largely deformed system, lul » 1 or lot » 1, first, the four oscillation modes are decoupled 
into two groups, two oscillations modes in the longitudinal directions, and the other two modes in the 
transverse directions. We can obtain these decoupled modes by taking C\ 2 = C14 = C 23 = C34 = 0. 
In the symmetric system, ri\ — n 2l furthermore, the two kinds of modes are further decoupled into the 
in-phase {Xl(t).i = ^l(t),2) an d the out-of-phase (Xl(t).i — ~ ^L(T).2) modes; namely the four excited 
oscillations are decoupled into the in-phase transverse breathing (ITB), in-phase longitudinal breathing 
(ILB), out-of-phase transverse breathing (OTB) and out-of-phase longitudinal breathing (OLB) modes. 

The oscillation frequencies in the symmetric system become 
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where the superscripts 'in' and 'out'' indicate the in-phase and out-of-phase modes, respectively, and the 
subscripts 'T' and 'L' denote the oscillations in the transverse and longitudinal directions, respectively. 
Note that f 1 =f 2 ,X 1 = X 2 and so on. 

Next we discuss the oscillations in the near-spherical symmetric system, lol ~ wr ~ 1- In this system 
the oscillation modes are decoupled into the monopole and quadrupole oscillation modes. 

By taking the collective coordinates to be Xla — Xt,i — K , we can obtain the frequencies of the 
monopole oscillations with the following equation: 
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By taking the collective coordinates to be At,-, 
frequencies of the quadrupole oscillations with 



= Xi, similarly, we can also calculate the 
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In the symmetric system pf^ = p% \ these modes are also decoupled to the in-phase and the out-of- 



phase modes. Then the oscillation frequencies in the symmetric system become 
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where the subscript 'M' and 'Q' indicate the monopole and quadrupole frequencies, respectively. 



(27) 
(28) 



III. RESULTS 

In this section we discuss the collective oscillations coupled with the transverse and longitudinal breath- 
ing modes. As mentioned in the previous section we fix the fermi energy to be ef = 20/27, which makes 
the parametric and ferromagnetic spin phase transition at g = 1 when B = 0. 

Density profiles illustrating these three regimes are shown in Fig. ^ The panels show the densities as 
a function of distance x in a weak magnetic field (g 2 Ns = 1.0 x 10~ 4 ) and at three different interaction 
strengths: g = 0.5(a), 0.95(b) and 1.05(c). The solid and dashed lines represent the scaled density distri- 
bution of major and minor components of fermions. From top to bottom, the panels show paramagnetic 
(a,b) and ferromagnetic (c). 

Next we calculate the collective oscillations by solving eigenvalue equation (121)1) . In this equation we 
obtain four eigenvalues corresponding to collective oscillation frequencies. For convenience we define these 
four frequencies as oj\, lu 2 , lo^ and W4 in order of the frequencies from the larger one to the lower one. In 
addition we refer to the four collective modes with the frequencies to 2 , U3 and W4 as modes-1, -2, -3 
and -4, respectively. 

In Fig. we show the frequencies of oscillations as functions of the deformation parameter cItl = 
]h(u)l/u>t) at the coupling constant g = 0.3 (a), = 0.6 (b) and = 0.9 (c). The solid lines represent the 
first and third frequencies io\ and W3, and the chain-dotted lines indicate lo 2 and U4. In these calculation 
we take the magnetic field to be g 2 B = 1.0 x 10 -4 , where the system is almost symmetric between the 
two component, N% ~ N 2 . In the same figures, furthermore, we also plot the frequencies of the transverse 
and longitudinal breathing modes with the dashed and dotted lines, respectively. 

First we see the strong level mixing between modes-2 and -3 at two deformation parameters, cItl = d\ 
and d 2 (d\ < < d 2 ). When drh <d\ and d 2 < dxL, the frequencies in the full calculation (solid lines) 
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are almost the same as those in the calculations decoupled with the transverse and longitudinal modes 
(dashed and dotted lines), though they are quite different when d\ < cItl < ^2- 

In Fig. we also give the same quantities but with g 2 B = 1.0 x I0~ 2 . The system is not symmetric 
Til 7^ n-2, but the results are quite similar to those in Fig. [2] though the difference in the frequency at the 
level mixing points between the modes-2 and -3 is a little larger than those in Fig. [21 

For future discussions we refer to the three deformation parameter regions, dxL < d\, d\ < dxh < <^2 
and c?2 < g?tl as the prolate deformed, the semi-spherical and the prolate deformed regions, respectively. 
The transverse and longitudinal breathing oscillation modes are almost decoupled in both the oblate and 
prolate deformed regions, while these modes are coupled in the spherical region. 

In order to know the level mixing behaviors more, we calculate the four following components: 

X't = ^(At,i + A Tj2 ), XT = "^( A i,i + (29) 

At"* = -^(At,i - A T , 2 ), Ar = j=(\l,i - Al, 2 ), (30) 

where Xp, Xf, X^ ut and X° L ut are components of the ITB, ILB, OTB and OLB modes, respectively. In 
Fig. 0] we show the components of the eigenvectors at the mode-1 (a), the mode-2 (b), the mode-3 (c) 
and the mode-4 (d) with g = 0.9 and g 2 B = 1.0 x 10~ 4 . The solid, chain-dotted, dashed and dotted lines 
represent A™, A™, X™ 1 and A£ ut , respectively. 

The results show that the phase of the modes-1 and -4 are in-phase in all drL regime. When the 
system is oblate deformed (g?tl << 0). the mode-1 is almost the ITB mode. As dxL increase, the system 
approaches the spherical, the ILB component in the mode-1 mode gradually becomes larger. When the 
system becomes spherical (g?tl = 0), the mode-1 becomes a mixed mode with the ITB and ILB in the 
same ratio: namely the mode-1 is the monopole oscillation. As the system is further deformed to the 
prolate shape, the transverse component decreases and the longitudinal one increases and the oscillation 
becomes the in-phase longitudinal oscillation in uj^/lot — ► oo limit. 

The mode-4 has the same behavior, but its phase is out-of-phase; the OTB mode when g?xl < rfi, the 
out-of-phase monopole mode at dxL — and the OLB mode when d-ri > <^2- 

The variations of the modes-2 and -3 show more complicated behaviors. The mode-2 is mainly the 
OTB mode in the oblate deformed region (c?tl < di), and the OLB mode in the prolate deformation 
region ( drL > ^2)- The mode-3 is the ILB mode in the oblate deformation, and the ITB mode in 
the prolate deformation. At dxL = di and d2, the modes-2 and -3 make the level mixing and exchange 
their roles. Namely the modes-2 and -3 are in-phase out-of-phase oscillations in the semi-spherical region, 
respectively. As dxL approaches zero, the TB and LB oscillation modes are mixed. At dxL = 0, the mode- 
2 becomes the in-phase quadrupole oscillations, (A™ = — 2A™) and the mode-3 becomes the out-of-phase 
ones, respectively. 

Thus we know that the four modes agree with the frequencies of the TB and LB modes, when lol/ojt >> 
1 (drL >> 0), and lut/^l >> 1 (drL << 0). On the other hand these modes show different behaviors 
in the nearly spherical symmetric system (drL ~ 0), where the collective oscillations are described as the 
coupled modes of the monopole and quadrupole oscillations. 

In Fig. El furthermore, we show the same quantities in Fig.0]but with g 2 B — 1.0 x 10~ 2 . The feature 
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discussed above is less clear, but also seen in the case of the strong magnetic field. 

When g 2 B = 1.0 x 10~ 2 , the number asymmetry is not so small, Ni > N2. The mixing between the 
in-phase and out-of-phase modes is very small in modes-1 and -4. The similar result was shown for the 
monopole oscillation in the spherical trap. Even in modes-2 and -3, this mixing is seen only in the region 
near cLtl ~ di and cLtl ~ d<i ■ Hence we can disregard the mixing between the in-phase and out-of-phase 
modes in qualitative discussions, and consider that modes-2 and -3 exchange their roles around cItl ~ di 
and dxL ~ di- 

In order to know more details of the deformation dependence, we should examine relations between 
the frequencies of the collective oscillations and the coupling constant g. First we show the frequencies 
of the monopole and quadrupole oscillations versus the coupling constant g in the spherical symmetric 
system (lul — lut) with the g 2 B = 1.0 x 10~ 4 in Fig.[S] The solid and dashed lines represent the results of 
the monopole and quadrupole oscillations. In the completely spherical symmetric system the monopole 
and quadrupole states are decoupled, and there is no level mixing between these two modes. In the 
paramagnetic phase (g < 1.0), larger frequencies correspond to the in-phase oscillation modes both in the 
monopole and quadrupole modes. The monopole modes make the level mixing between the in-phase and 
out-of-phase modes in the ferromagnetic phase {g < 1-0), while the quadrupole modes do not show the 
level mixing. As the coupling becomes larger, the frequency of the oscillation modes except the in-phase 
monopole oscillation mode decreases monotonically until g ss 1.0, and increases above that, while the 
frequencies of the in-phase monopole oscillation are almost invalid. When the coupling constant increases 
further, the level mixing between the two modes of the monopole oscillation occurs. 

Next we study the non-symmetric system, where all the four modes can be mixed. In Fig. \7\ we show 
the frequencies of the mode-1, -2, -3 and -4 as the functions of the coupling constant with lul/lut = 1.5 
{d TL = 0.41) (a), uj L /uj T = 0.7 {d TL = -0.36) (b), cj l /ujt = 0.5 (d TL = -0.69) (c) and uj l /uj t = 0.3 
(drL = —1-2) (d). In these figures we see again that the levels between modes-2 and -3 are mixed at 
a certain g. When uol/uot = 0.7 (Fig. Q), for example, the frequencies of modes-1 and -2 are u> = 2lot, 
and those of modes-3 and -4 are ui — 2ul at g — 0. As the coupling becomes larger, the frequencies of 
modes-2 and -4 monotonously decreases while those of modes-1 and -3 is almost invalid, and then the 
the level mixing between modes-2 and -3 occurs at g « 0.37. 

As the deformation parameter cLtl increases, the difference between lot and lol becomes larger, and 
the level mixing occurs at larger coupling constant. As the deformation parameter further increases, 
the level mixing disappears when cLtl < d\ because the minimum value of loi is larger than U3, and 
the level mixing cannot occur. The 0J2 becomes minimum at the border between the paramagnetic and 
ferromagnetic phases 5 ~ 1, where the level mixing occurs when <Itl — d\. 

In Fig. [H]we show the components of the eigenvectors at the mode-1 (a), the mode-2 (b), the mode-3 
(c) and the mode-4 (d) with g 2 B — 1.0 x 10 -4 and uil/uit = 0.7. The mode-1 and the mode-4 are almost 
ITB and OLB modes in the weak coupling region, respectively. As the coupling, g, increases, the ILB 
and OTB components are gradually mixed in both the modes. In the case of the mode-1, this mixing is 
small, while in the mode-4 the OLB and OTB components are more largely mixed, particularly around 

In this figure, furthermore, we can also see the level mixing between modes-2 and -3. At g = the 
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mode-2 and the mode-3 are the ILB and OTB modes. As the coupling increases, the ITB and OLB 
components are more mixed in the mode-2 and the mode-3, respectively. At g ss 0.37 the level mixing 
occurs. The strength of the components is exchanged between the modes-2 and -3. 



IV. SUMMARY 



In this paper we study the collective breathing oscillations in the deformed traps. The oscillations 
contain the four kinds of the oscillation, ITB, OTB, ILB and OLB modes. When the shape of the 
traps are largely deformed, either oblate or prolate, these modes are decoupled. As this shape becomes 
close to being spherical, the transverse and longitudinal modes are coupled, and make the monopole 
and quadrupole oscillations. The borders between these coupling and decoupling are at cItl ~ d\ and 
cLtl ~ d,2, where the OTB and ILB modes make level crossing. 

In the deformed traps the frequencies of the ITB and ILB modes agree with twice of the trap frequencies 
and are not so different from those of the non-interacting system. The frequencies of the OTB and OLB 
modes show some coupling dependence and give information of the interaction, Their qualitative behaviors 
are however similar to those in the dipole oscillation [jj. 

In the semi-spherical traps, on the other hand, the collective oscillations exhibit more complicated be- 
havior because of the level mixing. Hence we can get significant information on many body system of fcrmi 
gas by investigate breathing oscillations with the various fermion-fermion couplings and deformations of 
the trap near being spherical (0.4 < lul/lut ;$ 1-6). 

In this work we do not discuss extremely deformed system (wt - > oo or lul — > oo). In such extreme 
conditions the system becomes a quasi-low dimension one |l4| . and the collective oscillation can be 
expected to show quite different behaviors from those in our work 
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FIG. 1: The scaled density distribution of two components fermion at g = 0.5 (a), g = 0.95 (b) and g — 1.05 (c) 
with with e/ = 20/27 and g 2 B = 1.0 x 10 -4 . The solid and dashed lines represents the distribution of the major 
and minor components, respectively. 
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FIG. 2: The frequencies of the collective oscillations versus u)l/<-ot at g = 0.3 (a), g — 0.6 (b) and g = 0.9 (c) 
with e/ = 20/27 and g 2 B — 1.0 x 10~ 4 . The solid lines represent the results with the full calculations, and the 
dashed and dotted lines indicate the frequencies of the transverse and longitudinal breathing modes, respectively. 
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FIG. 3: Same as Fig.H but with g 2 B = 1.0 x 10~ 2 . 
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FIG. 4: Components of the eigenvectors of mode-1 (a), mode-2 (b), mode-3 (c) and mode-4 (d) at g — 0.9 with 
g 2 B = 1.0 x 10~ 4 . The solid and chain-dotted lines represent the in-phase transverse and longitudinal components, 
respectively, and the dashed and dotted lines denote the out-of-phase transverse and longitudinal components, 
respectively. 
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FIG. 5: Same as Fig.g] but with g 2 B = 1.0 x 10" 




FIG. 6: The frequency of the monopole and quadrupole oscillations with the external magnetic fields with g 2 B = 
1.0 x 10 -4 . The solid and dashed lines represent the results of the monopole and quadrupoles modes, respectively. 
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FIG. 7: The frequencies of the collective oscillations versus the coupling constant g with the external magnetic 
fields with g 2 B = 1.0 x 10~ 4 for lu l /cj t = 1.5 (a), 0.7 (b), 0.5 (c) and 0.3 (b). The solid lines represent the results 
of the mode-1 and the mode-3, and the chain-dotted lines indicate those of the mode-2 and the mode-4. 
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FIG. 8: Components of the eigenvectors of mode-1 (a), mode- 2 (b), mode-3 (c) and mode-4 (d) with g 2 B = 
1.0 x 10~ 4 and lul/^t = 0.7. The solid and chain-dotted lines represent the in-phase transverse and longitudinal 
components, respectively, and the dashed and dotted lines denote the out-of-phase transverse and longitudinal 
components, respectively. 



